Mon. Not. R. Astron. Soc. 000, [HH] (2008) Printed 19 May 2009 (MN MrgX style file v2.2) 



First Results from the CHEPS^: Exoplanets and the 
Discovery of an Eccentric Brown Dwarf in the Desertf 

J.S. Jenkins^, H.R.A. Jones^, K. Gozdziewski^, C. Migaszewski^, J.R. Barnes^, 
M.I. Jones\ P. Rojo\ D.J. Pinfield^, A.C. Day-Jones^, S. Hoyer^ 

^Department of Astronomy, Universidad de Chile, Casilla 36-D, Santiago, Chile, email: jjenkins&das.uchile.cl 
, ^ Center for Astrophysics, University of Hertfordshire, College Lane Campus, Hatfield, Hertfordshire, UK, ALIO 9AB 
^Torun Centre for Astronomy, Nicolaus Copernicus University, Gagarina 11, 87-100 Torun, Poland 

• Draft: 11/08 



ABSTRACT 

We report the discovery of a brown dwarf on an eccentric orbit and with a semima- 
jor axis that places it in the brown dwarf desert region around the star HD191760. The 
star has a spectral type of G3IV/V and a metallicity ([Fe/H]) of 0.29 dex. HD191760 
adds to the small number of metal-rich stars with brown dwarf companions. The 
brown dwarf (HD1917605) is found to have an orbital period of 505.57±0.40 days 
and semimajor axis of 1.35±0.01 AU, placing it firmly in the brown dwarf desert. 
The eccentricity of HD191760& is found to be 0.63±0.01, meaning it reaches as close 
as 0.5 AU from the host star. Dynamical simulations indicate that no inner planets 
could reside at separations beyond ^0.17 AU due to the disastrous gravity imposed by 
HD191760&. In addition to these first results we also refine the orbits found for the ex- 
oplanets around the stars HD48265, HD143361 and HD154672. All 1-planet solutions 
are in agreement with those previously published by the Magellan Planet Search. 
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1 INTRODUCTION 

Since the first few extrasolar planets (hereafter exoplan- 
ets) were detected (iMavor fc _Queloj[l995l : lMarcv fc Butled 
Il996al : iMarcv fc Butleiill996bl ) it became readily apparent 
that planets preferred metal-rich environments (iGonzalej 
1 1997^ ) . This discovery has been placed in the framework of gi- 
ant planet formation through core accretion of gas depleted 
materials left over from the formation o f the parent star (e.g. 
iKornet et alll2005l : lAlibert et al.ll2005h . This hypothesis has 
been extensively tested by different methods and an alysis 
techniques, most recently bv lFischer fc Valentil |200^, and 
remains the most plausible scenario for the increased proba- 
bility of planet detection around metal-rich stars. The com- 
peting theory is based on the concept that planets migrate 
through the proto-planetary disks from wh ich they formed 
l|Lin fc Papaloizoull 19861 : iTrilhng et aLll2002l ') . It has been hy- 



* Calan-Hertfordshire Extrasolar Planet Search 
t Based on observations collected at the La Silla Paranal Obser- 
vatory, ESO (Chile) with the HARPS spectrograph on the ESO 
3.6m telescope, under the program IDs 079.C-0927(B), 079.C- 
0927(C), 081.C-0148(A), 081.C-0148(B) and 282.C-5034. 



pothesized that the stellar atmosphere is polluted after for- 
mation through i nfall of migrating pl a nets o nto the stellar 
surface. However, iPinsonneault et al.l l|200ll ) modelled the 
infall of material onto stellar atmospheres at the ZAMS and 
found that 10 Me worth of gas depleted material could raise 
the metallicities ([Fe/H]) of G dwarfs to ~0.3 dex, whereas 
the same material falling on an F dwarf would increase the 
metallicity to ~0.6 dex. In their large sample (>1000 stars) 
of consistent metallicities, Fischer fc Valenti found no cor- 
relation between the depth of convective zones of exoplanet 
hosts and their metallicity, indicating that infall of planets 
or gas depleted material can not explain the over abundance 
of metal-rich atmospheres for exoplanet host stars. 

Whatever the scenario that leads to exoplanet host stars 
being rich in metals, it is clear that such stars have a higher 
probability of hosting exoplanets. Fischer fc Valenti have de- 
termined that the fraction of planets increases with increas- 
ing metallicity proportional to the square of the number of 
iron atoms in the stellar photosphere. This correlation is be- 
ing exploited by a number of planet search projects like the 
N2k, Keck, etc who bias their samples, or subsets thereof, 
towards the most metal-rich stars. We are exploiting this 
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correlation also, by searching around the most metal-rich 
solar-type stars in the southern hemisphere as part of our 
CHEPS projec t The CHEPS targe t list is compiled from 
the analysis of Ijenkins et al] l|2008l ') by utilising high S/N 
high resol ution spectra, acq uired using the FEROS spec- 
trograph (jKaufer et al.l[l999h . The targets were selected to 
have [Fe/H] ^0.1 dex and all have B — V colours in the range 
0.5-0.9 (Spectral Types of late-F to mid-K). Along with the 
metallicity bias we are also focusing on the most inactive 
stars. The Ca ll HK line cores were used in Jenkins e t al. to 
measure the logiinK activity index l|Noves et al ]|1984| ) which 
has been shown to be a useful proxy to the level of noise 
(aka jitter) in radial-velocity times eries (ISaar fc Donahu3 
119971 : isaar eTaLlllQQSl : IWrighdl2005h . We therefore selected 
our primary sample to include all stars with logi?HK^"4.5, 
which allows for the lowest levels of radial- velocity jitter in 
our dataset. 



2 OBSERVATIONS 

We present data fro m six observing ru ns using the ESO- 
HARPS instrument (|Mavor et al.1 120031 ) on the 3.6m tele- 
scope in La Silla Chile. HARPS is a fiber-fed cross-dispersed 
echelle spectrograph. The instrument employs two 1" on sky 
fibers, one on object and one on the Thorium- Argon (ThAr) 
calibration source, which feed the light to the spectrograph 
to be re-imaged on two 2k4 CCD chips. A total of 72 orders 
are spread over the two CCD's, 36 object and 36 calibration 
source, with an instrumental resolution [R) of ~115,000. 29 
observations were made of HD191760, with 8 for HD48265 
and 6 each for both HD143361 and HD154672, all using 
our standard observational procedure, which aims to ob- 
serve each star for a fixed integration time of 15 minutes 
in order to integrate over the strongest p-modc oscillations 
that all solar-hke stars exhibit (lO'TooTe et al.ll2008l ). Also to 
attempt to increase our efhciency to the detection of short 
period exoplanets, and in particular low mass exoplanets, we 
intensively monitor the stars on short timescales to search 
for potential planetary signals, including for HD191760. 

The reduction and analysis of all stars are performed 
in near real time at the telescope, us ing the latest version 
of the HARPS-DRS (|Pepe et al.ll20oi ). The DRS (Data Re- 
duction Software) performs all functions of the reduction 
and analysis, from bias removal, order localisation, flatfield- 
ing, cosmic ray removal, scattered light subtraction, extrac- 
tion and blaze correction, to determination of the radial- 
velocities, barycentric Julian dates and internal uncertain- 
ties through cross-correlation with stellar templates and the 
reference ThAr calibration orders. Due to the high quality 
mechanical, thermal and pressure stabilisation of HARPS, 
the drift over a single night is found to be less than 1 m/s. 
Such a system provides internal uncertainties below 1 m/s, 
meaning our short term p recision is domin ated by stellar 
activity induced jitter (e.g. ISaar et al.iri998l ). 



3 HD191760 CHARACTERISTICS 



(|Perrvman et al ] | 19971 : Ivan Leeuwenll2007h and also with the 
derived distance fro m the Hipparcos main sequence (AM^; 
'Jen kins et al.|[200^ . All Hipparcos astrometry was obtained 
using the AstroGrid tool which lists the new catalogue val- 
ues from van Leeuwen with i ncreased astrometric pre cision. 
The kinem atic catalogue of 
updated bv lHolmberg et ahl 



Nordstrom et al.1 (|2004l '). later 
20071 ). includes observations of 



Table [T ]lists the stellar characteristics of the star HD191760. 
iHoukI (|l978l ') lists the star's spectral type as G3IV/V, 
which is in agreement with the value adopted by Hipparcos 



HD191760 through Stromgren filters and ca librations based 
on the infrared fiux method (|Alonso et al.iri996i l to derive 
various stellar characteristics. They also used parallax mea- 
surements from Hipparcos, alo ng with Padova evolution- 
ary tracks l|Girardi et al.l |2000| ) to determine evolutionary 
parameters for their sample. We have also derived vari- 
ous parameters for a number of stars in preparation for 
our metallicity based planet search, both by high resolu- 
tion FEROS spectra and precise photometry. We used the 
Hipparcos optical magnitudes, along with 2MASS photome- 
try, to measure the eff ective temperatures of our stars by the 
infrared flux method (jBlackwell et al.lll990l ). We used these 
temperatures to interpol ate a grid of high resolution Ku- 
rucz model atmospheres l|Kurucj|l993l ). built using WITA6 
(Pavlcnko 2000), to de termine the stellar met allicity. Yonsei- 
Yale (Y2) isochrones (iDemargue et al]|2004l) were then in- 
terpolated using these temperatures, metallicities and an es- 
timate for the alpha-enhance ment parameter that wa s ob- 
tained by the locus of points in lEdvardsson et al] (|l993l ) and 
iPagel fc Tautvaisiend (|l995l ). 



The properties of HD191760 derived by Holmberg et 
al. and this work mostly agree within the formal uncertain- 
ties. Their metallicity value of 0.12±0.10 dex is fou nd to be 
smalle r than our value of 0.29±0.07 dex and .Holmberg et al.l 
(|2008l ) show that there is a small systematic offset be- 
tween their metallicities and those published in Jenkins 
et al. due to the differing effective temperature scales em- 
ployed. Fig. [T] shows the position of HD191760 on an HR- 
diagram, along with the Y2 isochrones in steps of Q.IMq 
(solid curves) and higher resolution steps of O.OIM© (dashed 
curves). The position of the star, given the bolometric lu- 
minosity of 2.69±O.2L0, effective temperature of 5821±82K 
and [Fe/H] of 0.29±0.07 dex, is best explained by a star of 
mass 1.281q'5o-^0, where the uncertainties are dominated 
by the uncertainty on the Hipparcos parallax. Our estimate 
for the age of the star is 4.ll2'gGyrs. The large lower er- 
ror bars arise due to the position of the star close to the 
main sequence turnoff. The surface gravity (loggr) value of 
4.131q Q4dex is consistent with a star evolving onto the sub- 
giant branch. 



As part of the analysis we performed to select our tar- 
get sample we also extracted chromospheric activity indices 
from measurements of the level of emission in the Ca ll HK 
lines. The logi?'HK-index has been shown to correlate well 
with the leve l of noise (a ka jitter) in precise radial-velocity 
datasets fe.g. lSantos et al . 2000; Wright 2005). Fig. [2] shows 
the Ca II H line core of HD191760 (top) along with the other 
three planet-host stars we present data for in §5. Our low 
derived value of log7?'HK=-5.17 is highlighted from the deep 
central line core, or lack of core re-emission, of HD191760. 
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Figure 1. Evolutionary Y2 mass tracks are shown in steps of 
O.IM0 (solid curves) and finer steps of O.OIMq (dashed curves) 
for effective temperatures and bolometric luminosities. HD191760 
is represented by the error bars and its position is consistent with 
a mass of 1.28lo'5oM0 and an age of 4.21I2 s^yrs. The masses 
of the solid curves are shown for reference. 



Table 1. Stellar parameters for HD191760. 



Parameter 



HD191760 



Spectral Type 

logfl'HK 
Hipparcos N^^^ 
Hipparcos cr 
TT (mas) 

M^/Mq - H07 
M»/M0 - Here 

R*/Rq 
T,s (K) - H07 
Toff (K) - JOS 
[Fe/H] - H07 
[Fe/H] - JOS 

log(9) 
U,V,W - H07 (km/s) 
vsinij^i (km/s) 
tisiniHOT (km/s) 
DsinicCF (km/s) 
Prot R' (days) 
^Vot.CCF (days) 
Age^j^^ (Gyrs) 
Age (Gyrs) - H07 
Age (Gyrs) - Here 
Jitter - logi?'HK (™/s) 



G3IV/V 

-5.17 

99 

0.011 

12.25±0.90 

1.107 

2.69±0.20 
1 14+0.08 

1 90+0.02 

1.62±0.07 
5794±76 
5S21±82 
0.12±0.10 
0.29±0.07 

4 1Q+0.05 
^•^^•^-0.04 

-29,-21,15 
3.47 
2 

2.33±0.05 
25.2 
35.1 
9.9 
5 fi+O-9 
4 1+0-8 

^•^-2.8 
1.50 
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V) 

E, 
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H07 relates t olHolmberg et al. JOS is the reference for 

I Jenkins et al ] l|2008h and Here labels values derived in this work. 



4 ORBITAL SOLUTION FOR HD191760 

The 29 Dopper velocities for HD191760 were measured in 
six observing runs, four of which spanned five nights in du- 
ration, with the final two points acquired through ESQ Di- 
rectors Discretionary Time. All velocities and bisector span 
values (BVS) are listed in Tabled note that the two large 
deviant BVS values are from our two lowest quality spec- 
tra which were interrupted by clouds, significantly lowering 
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Figure 2. The calcium o H line cores for the four stars discussed 
in this work. From top to bottom we show HD191760, HD4S265, 
HD143361 and HD154672 respectively. All four objects exhibit 
no re-emission of flux in the line core, highlighting the inactive 
nature of these stars. 
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Figure 3. Panel (a) shows the best fit Keplerian solution to 
the Dopper points fo r HD 191760 with BVS correction as in 
iMigaszewski fc NowaQ ll200a) . The data come from six observing 
runs spread over a full orbital period of the companion. The addi- 
tion of two Doppler points close to periastron passage was used to 
confirm the eccentric nature of the companion. Panel (b) shows 
the residuals without BVS correction and panel (c) shows the 
residuals after correction, corresponding to the synthetic curve in 
panel (a). Panel (d) shows the radial- velocities against the BVS 
after subtraction of the best fit Keplerian synthetic curve. Panel 
(e) shows the best fit linear correlation coefficient of BVS-RV. 
Note that the two values with the largest BVS help to better 
constrain the correlation. 



the S/N, meaning the construction of the CCF was much 
poorer than the rest of the high S/N spectra acquired. Panel 
(a) in Fig. |3] shows these velocities with the best fit Ke- 
plerian solution to the data. Due to the scheduling of the 
runs the orbit is firmly constrained by the final two data 
points. Before the final two data points were acquired, a 
resonant double planet solution could describe the dataset 
just as well. Such system ambiguities could be common 
in radial-velocity datasets, p articularly at low eccentricities 
(lAnglada-Escude et ahllioOSl ). However, in this case the sys- 
tem configurations were so contrasting that one data point 
was enough to decouple them, and with an eccentricity of 
around 0.63 and no indication of any trend in the third ob- 
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Table 2. Radial-velocities for HD191760 



BJD 


RV 


Uncertainty 


BVS 


RV-aBVb 


(-2454000) 


(m/s) 


(m/s) 


(m/s) 


(m/s) 


250.90047 


-30662.99 


0.88 


9.49 


-235.23 


253.89490 


-30689.15 


0.70 


3.20 


-257.01 


365.50781 


-30575.97 


1.37 


9.45 


-148.19 


365.73781 


-30577.79 


0.62 


5.68 


-147.39 


367.49679 


-30540.86 


0.38 


3.66 


-109.04 


367.72155 


-30535.38 


0.73 


-0.09 


-100.95 


368.49981 


-30519.40 


0.53 


-0.94 


-84.39 


368.71411 


-30498.29 


1.86 


30.70 


-85.31 


369.49468 


-30503.67 


0.77 


-1.71 


-68.12 


369.62706 


-30501.89 


0.50 


0.91 


-68.16 


369.68010 


-30497.67 


0.51 


2.50 


-65.07 


369.74838 


-30497.94 


0.63 


-0.93 


-62.93 


577.77658 


-30045.01 


0.82 


7.41 


384.19 


577.90446 


-30051.42 


0.40 


5.55 


379.07 


578.77151 


-30047.07 


0.54 


6.49 


382.77 


578.90713 


-30049.68 


0.48 


6.46 


380.18 


579.74685 


-30048.12 


0.56 


3.12 


384.07 


579.90348 


-30049.27 


0.41 


5.07 


381.57 


580.74769 


-30052.14 


0.60 


6.86 


377.44 


580.89930 


-30051.97 


0.51 


5.94 


378.25 


581.76210 


-30049.85 


0.39 


5.45 


380.71 


581.91331 


-30052.58 


0.36 


7.26 


376.73 


723.53335 


-30373.95 


3.17 


93.78 


-4.90 


724.55839 


-30440.36 


0.34 


7.46 


-11.19 


725.53823 


-30446.88 


0.44 


8.10 


-18.15 


726.53101 


-30452.77 


0.43 


7.76 


-23.81 


727.53134 


-30458.24 


1.09 


1.25 


-24.74 


812.52135 


-31646.47 


1.02 


12.07 


-1220.52 


813.53608 


-31679.70 


1.16 


10.27 


-1252.50 



servational dataset, which is where a resonant system should 
most manifest a trend, we believe the quoted solution is ac- 
curate. Panel (b) shows the residuals to this fit without cor- 
recting for the bisector variations. Such variations have been 
shown to be a robust tracer of line asymmetries that can 
impa c t stellar radial-ve locity measurements ( Queloz et al.l 
I2OOII : iHenrv et al.ll2002l ). We perfor m corrections for such 
variations following the analysis in iMigaszewski fc NowakI 
l|2009l ) and the residuals for the data after BVS correction 
is shown in panel (c). The radial- velocities (RV) are shown 
against the bisector velocity span (BVS) values in panel (d). 
This hightlights the difference of the two offset velocities, 
particularly one data point which has a BVS value greater 
than 90m/s. The correlation coefficient (a) used to correct 
the radial-velocities is shown in panel (e) by the solid line. 
It is clear that the two velocities with large BVS values help 
to constrain the correlation across the parameter space and 
indicate the technique is robust, even for fairly large BVS 
offsets. The a correction factor is found to be 0.697±0.064 
and the corrected radial- velocities, those that are plotted in 
panel (a), are listed in column 5 of Table. [2] Clearly the data 
describes the fit much more accurately after BVS correction 
has been performed, compared with no corrections shown in 
panel (b). Given the level of the residuals, with uncertainties 
which include in quadrature the stellar jitter of 1.50 m/s, be- 
ing all around ilcr from the fit, the quoted orbital solution 
is accurately constrained. The orbital characteristics for a 
single companion solution is shown in column 2 of Table, [l] 
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Figure 4. Stability maps for a potential planet in semi- 
major axis— eccentricity space in the HD191760 system. A pos- 
sible 2-planet solution of reasonable quality and low eccen- 
tricity of the inner companion was found, and is shown by 
the cross. The top panel is for the S pectral Number indica- 
tor llMichtchenko fc Ferraz-Mellol I2OOII ) , the lower panel is for 
the maximal eccentricity attained by the putative planet dur- 
ing 10^ yrs of integration time. The scale at the top of the panels 
represents the zones of stability, with the darkest regions the most 
stable. Both indicators around the possible planetary object show 
no zones of stability anywhere. A sharp border of stable motions 
is visible. Any planetary companion is unlikely to exist beyond 
0.17-0.18 AU regardless of its initial eccentricity. 



5 POTENTIAL FOR FURTHER COMPANION 
OBJECTS 

In order to map out a strategy for future radial-velocity ob- 
servations of this star, or indeed to decide if further study is 
warranted, we decided to look for additional hidden compan- 
ions in the data. This seems plausable given that the signal 
describes a fairly wide orbiting companion near the plan- 
etary mass boundary, therefore it is reasonable to assume 
there are other companions in the system with significantly 
lower mass that would induce a secondary orbital signature 
with a lower amplitude. Therefore we performed a search 
for additional companion planets in the system even though 
the current dataset is rather limited. 

First, we tested a hypothesis that the radial- velocity sig- 
nal may be produced by an additional planet with an inter- 
mediate period orbit. We performed the search for additional 
planets using the genetic algorithm with MEGNO penalty 
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(GAM P) (jGozdziewski et aljlioo^ : iGozdziewski fc Konackil 
I2OO6I ). Essentially, this numerical approach searches the 
phase space using a genetic algorithm to find zones of sta- 
bility where additional companions could reside. Unfortu- 
nately, the quite extensive GAMP search did not bring 
any reasonable stable solutions beyond '^0.2 AU limit. Ac- 
tually, we might expect that relatively far from the star, 
only strongly resonant configurations could be found due to 
strong gravitational influence of the brown dwarf. We tested 
this by computing dynamical stability maps for reasonable 
fits found without considering any stability checks. Such an 
example configuration returns a secondary planetary object 
with a minimum mass of 0.10 Mj, a low eccentricity of ~0.17 
and a semimajor axis of ^^0.22 AU. Such a configuration 
gives rise to a reduced rms of 1.84 m/s and could be hidden 
in the current dataset. However, as shown in Fig.|4]the entire 
region around the hypothetical planet is extremely chaotic 
and therefore it is extremely unlikely that such a planet 
could reside at this location. The two panels of this figure 
represent stability maps in the semimajor axis-eccentricity 
plane in terms of the so called Spectral Number indicator 
of fo rmally chaotic motions (Michtch^nko & Fcrraz-Mello 
I2OOII . top panel) and the maximum eccentricity indicator 
(the right-hand panel), i.e. the maximum eccentricity at- 
tained during the integration period indicating the degree of 
geometric instability (large max Cc means collisions). Both 
maps were computed over 10^ years that corresponds to 
^ 50, 000 revolutions of the outer companion. The figure 
shows that we must rule out, not only this potential compan- 
ion, but also all companions in the vicinity of this parameter 
space. Actually, the strongly unstable zone ends at ~ 0.17- 
0.18 AU, and it is unlikely any other planet could survive 
the disastrous, short-term perturbations of the brown dwarf 
companion beyond that limit. 



6 DISCUSSION 

6.1 Brown Dwarf Companion 

The constrained orbital parameters for HD1917606 place 
the minimum mass of the companion in the brown dwarf 
regime (M sin j=38.17±1.02Mj). Also the orbital semima- 
jor axis is found to be 1.35 AU, well inside the brown dwarf 
desert region. This so called brown dwarf desert was coined 
due to the paucity of such objects around solar-like star s 
with semimajor axes of <3-4 AU (|Marcv fc Butleill2000h . 
They showed that fewer than 0.5% of surveyed stars within 
3 AU on Doppler programs at the time had brown dwarf 
companions, whereas 5% of these stars harbour exoplanets. 
Therefore, since our sample is small (100 stars) we may not 
statistically expect to detect any brown dwarf desert com- 
panions and this object may have an actual mass above the 
hydrogen burning limit (~80 Mj) after inclination effects 
are considered. For stars with randomly distributed incli- 
nation s, 68% are found to fall within 47° (|Marcv &: ButleJ 
Il996al ). giving rise to a true mass of M=52.19Mj, comfort- 
ably below the hydrogen burning limit. An inclination i of, 
or below, 28.5° is required in order for HD191760fe to be a 
true stellar companion and such a scenario has an a priori 
probability of only 12.1%. Also there is no bin arity flag as of 
the latest update of the Hipparcos astrometry (|van LeeuwenI 



I2OO7I ). In addition, a visual inspection of the HARPS spec- 
tra does not reveal any evidence for additional stellar lines 
from an unseen companion and there was no extension of the 
acquisition image at the telescope to indicate another stel- 
lar component, ruling out the possibility of extremely high 
inclinati ons. However, such a scenario wou ld agree with the 
results of lStamatellos fc WhitworthI (|2009l ) who suggest that 
brown dwarfs form by disk fragmentation in such systems 
and that short period radial-velocity companions are prefer- 
entially low mass stars. This scenario would agree well with 
the high orbital eccentricity of the syst em and could also 
explai n a recent similar discovery from IWittenmver et all 
(2009). Interferometry could be used to help rule out a lower 
mass stellar companion. 

Another way to probe the inclination is to study the 
rotation velocity along with the activity index from the 
CaHK lines. The rotation velocity of this star is found to be 
2. 33±0.05 km/s , usin g a similar method to that explained 
in lSantos et al.l l|2002l ) (the uncertainty is the standard de- 
viation of the 25 measured cross correlation full widths at 
half maximum), which for the determined radius we find 
gives an upper limit to the rotation period of 35 days. The 
estimated r otation period from the log-RjjK-rotation period 
correlation (|Noves et al.lll983 ) is 25.2 days, which is incli- 
nation independent and should represent the equatorial ro- 
tation period, assuming that the orbital and rotation axes 
are aligned in the same plane. Assuming the rotation period 
from the activity relation holds then we estimate the system 
has a sin i=0.67, which gives an inclination of ~42°. Such 
an inclination would give rise to an actual mass of 57 Mj, 
which is still comfortably below the hydrogen burning limit. 
However, the log_RHK activity index was formulated for stars 
on the main sequence and may be gravity dependent so 
would not be fully applicable to this star. Also the use of 
this index should be made only after a number of stellar cy- 
cles have been observed, which allows one a robust estimate 
of the mean log-RR K- We have applied the value quoted in 
I Jenkins et al.|[2008l (one measurement), not the mean, so we 
do n ot know th e spread, which can vary at the ±0.1 dex level 
(see ljenkins et al. 2006), and can not apply this to the un- 
certainty on the relationship itself, which Noyes et al. quote 
as 0.08 dex for the convective turnover time. 

6.2 Benchmarkability? 

Since HD191760 is found to be evolving off the main se- 
quence onto the sub-giant branch the age estimates f or the 
system can be fairly well constrained. Holm berg et al.l (|2007l ) 
estimate an age of 5.6l5 g Gyrs in good agreement with our 
estimated value of 4.ll2'g Gyrs, the main differences arising 
due to the differing metallicities, and hence masses, found for 
this star. The mean of the two estimates gives rise to an age 
of 4.9±1.1 Gyrs, clearly a fairly old and evolved star. Such 
an age means that this system could provide a benchmark, 
evolved system, to test brown dwarf evolutionary models. 
Given there is a lack of well constrained age estimates for 
the older population of brown dwarf stars, if one can at- 
tain photometry or spectroscopy of the brown dwarf in this 
system then this could provide a well constrained metallicity 
and age to use in compari son with current evolutionary mod- 
els of brown dwarfs (e.g. iBaraffe et al.l [20031 : iBurrows et al.l 
I2OO6I ). Again interferometry may allow such studies to be 
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Table 3. Radial-velocities for HD48265, HD143361 and 
HD154672 



BJD (-2454000) 


RV (m/s) 


Uncertainty (m/s) 




HD48265 




365.824572 


23394.21 


3.47 


366.874164 


23399.40 


3.50 


367.852214 


23395.69 


3.45 


580.522907 


23449.18 


3.45 


581.553905 


23450.06 


3.45 


724.833698 


23413.42 


3.45 


725.817765 


23414.08 


3.44 


726.779101 


23411.87 


3.47 




HD143361 




250.745318 


-437.90 


2.39 


253.769447 


-436.44 


3.23 


367.556628 


-474.74 


2.22 


368.537299 


-474.20 


2.24 


578.788567 


-534.13 


2.18 


581.807762 


-531.75 


2.19 




HD154672 




248.713853 


-3796.59 


4.79 


250.731312 


-3794.42 


2.25 


367.595129 


-3653.29 


2.15 


368.618599 


-3651.00 


2.22 


578.811238 


-3797.81 


2.14 


581.870548 


-3808.85 


2.13 



performed on this system. Once the inclination is pinned 
down and the nature of the companion is confirmed as sub- 
stellar then we can also use the HARPS spectra to obtain a 
furt her age constrain t from the strength of the lithium line 
(see lMelo et al.|[2006l ). This would provide an additional, and 
independently determined, age estimate to further constrain 
the system. In addition, further atomic abundances could be 
garnered for the star to aid in constraining the metallicity 
and abundance nature of brown dwarf model atmospheres. 
Such an abundance study is already underway. 



7 REFINED ORBITS FOR THE PLANETS 
AROUND HD48265, HD143361 AND 
HD154672 

Due to some overlap in our target sample with that of 
the Magellan Planet Search, we can update a few of the 
recent planet discoveries announced by this program with 
additional velocities we have found at HARPS. Table. [3] 
shows our absolute Dopper velocities for the three planet 
host stars HD48265, HD143361 and HD154672. For aU of 
these stars we have acquired observations in three observ- 
ing runs with overall baselines around one year in duration. 
Fig. [3 shows the best fits to all three stars using Systemic 
l|Maschiari et al.ll2009l ), with the Magellan velocities in black 
(blue) and our HARPS values in grey (greenfl 

We show the updated orbital solutions in columns 
3, 4 and 5 of Table. |4] and in ge neral they agr e e wel l 
with the published results in both iMinniti et al.l (|2009| ) 



Colours in brackets are for the online version of the paper 




2453746.9 1 2454025.51 



Figure 5. The radial-velocity Keplerian fits to the stars HD48265 
(top), HD143361 (middle) and HD154672 (bottom). Plotted in 
blue are the literature values from the Magellan program and in 
green are the new data points from here. 



(HD48265 and HD143361) and lLopez-Morales et all (|2008l ) 
(HD154672). Lopez-Morales et al. used 16 velocity points in 
their orbital solution for HD 154672 over a number of orbital 
periods of the planet. We find a very similar solution with 
a slightly reduced rms of 3.37 m/s. Note that we used the 
stellar masses quoted by this group for a direct comparison 
of the results found. 

For HD48265 we find fair agreement with the derived 
properties assuming a 1-planet Keplerian solution. When 
adding our eight doppler velocities to the 17 acquired at 
Magellan, we find that the orbital period is now 700±8 days, 
albeit with a slightly increased rms of 6.78 m/s. The xl also 
increases when compared to the Magellan data alone from a 
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value of 1.37 up to 3.33. Given the baseline between the two 
data sets this might argue for an additional companion in the 
system, however currently both a search in the periodogram 
of the residuals to the 1-planet fit and the addition of a linear 
trend do not reveal anything significant. Indeed, although 
this appears to be the best fit solution, there is no significant 
power anywhere in the periodogram when combining the two 
data sets. 

Finally, for HDf 433616 we find good agreement with the 
published results, again with a slightly reduced orbital pe- 
riod of 1057 days and eccentricity of 0.15. This relates to an 
M sin i of 3.12 Mj, in good agreement with the 3 Mj found 
by Minniti et al. and with a lower rms of 3.37 m/s. For a 
1-planet solution none of these three systems exhibit signif- 
icant trends in the residuals and hence present no evidence, 
as yet, for additional longer period planets in the system. 
Note the uncertainties shown in Tab le. [3] are derived from 
the updated relation of I Wright l|2005h (private communica- 
tion) and the low level of each is highlighted in Fig. [2] where 
no Ca II H line core emission is seen in any of these three 
stars (HD48265: top-middle, HD143361: lower-middle and 
HD154672: bottom). 



8 SUMMARY 

We announce the first four radial-velocity detections from 
the CHEPS on the ESO-HARPS instrument. The high- 
light is the discovery of an eccentric brown dwarf in the 
desert orbiting the metal-rich G3IV-V star HD191760. For 
the best fit solution, HD1917606 is found to have a mass of 
38.17±1.02 Mj, an orbital period of 505.65±0.42 days and 
an eccentricity of 0.63±0.01, and this discovery adds to the 
small list of metal-rich stellar-brown dwarf binary systems. 
Our dynamical simulations show that no additional planets 
could exist inside the orbit of the brown dwarf with semima- 
jor axes larger than around 0.17 AU due to the destructive 
tidal forces of the brown dwarf. 

Finally, we also confirm the exoplanet detections of 
HD482656, HD1433616 and HD1546726 from the Magel- 
lan program by employing 1-planet Keplerian solutions to 
the data. We show updated orbits by adding eight veloc- 
ities for HD48265 and six velocities each for HD143361 
and HD154672. All orbital param eters are in agree - 
ment with the published results by iMinniti et al.l (|2009| ) 
and iLopez-Morales et all (|2008l '). Both HD1433616 and 
HD154672b show lower rms by the addition of our data, how- 
ever the best fit for HD48265ti shows an increased rms and 
xS. Since the star is shown to be chromospherically inactive 
this indicates that additional radial-velocities are required 
to probe for further companions in this system. 
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